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ABSTRACT 

Two groups have recently discovered a statistically significant deviation in the fluxes of high-redshift 
type la supernovae from the predictions of a Friedmann model with zero cosmological constant. In this 
Letter, I argue that bright, dusty, starburst galaxies would preferentially eject a dust component with 
a shallower opacity curve (hence less reddening) and a higher opacity/mass than the observed galactic 
dust which is left behind. Such dust could cause the falloff in flux at high-z without violating constraints 
on reddening or metallicity. The specific model presented is of needle-like dust, which is expected from 
the theory of crystal growth and has been detected in samples of interstellar dust. Carbon needles with 
conservative properties can supply the necessary opacity, and would very likely be ejected from galaxies 
as required. The model is not subject to the arguments given in the literature against grey dust, but 
may be constrained by future data from supernova searches done at higher redshift, in clusters, or over 
a larger frequency range. 

Subject headings: cosmology: observations - dust, extinction - radiative transfer 



1. INTRODUCTION 

Recently, two separate groups have interpreted their ob- 
servations of type la supernovae as evidence for acceler- 
ation in the cosmic expansion, presumably caused by a 
nonzero cosmological constant (Riess et al. 1998, hence- 
forth 'R98', and Perlmutter et al. 1998). Using the super- 
novae as standard candles (after correction for the relation 
between luminosity and light-curve shape), both groups 
find a progressive dimming of supernovae at high red- 
shift relative to the predictions of a flat, matter-dominated 
model, or even an open model with zero cosmological con- 
stant. There are at least two explanations for this dimming 
which are unrelated to the cosmological parameters: evo- 
lutionary effects and dust. The two observational groups 
have expended considerable effort attempting to account 
for such systematics, but this letter argues that the ob- 
scuration of distant supernovae by a cosmological distri- 
bution of dust is not ruled out, by developing a simple 
intergalactic dust model which is reasonable in its cre- 
ation and dispersement, has reasonable physical proper- 
ties, and could cause the observed effect without violating 
constraints such as those noted in R98. 

The standard way to estimate dust extinction, used by 
the supernova groups, is to apply a relation between red- 
dening and extinction derived from Galactic data (see, e.g., 
Cardelli, Clayton & Mathis 1989). This relation is a result 
of the frequency dependence of the opacity curve of the ab- 
sorbing dust, which can be measured well for Galactic dust 
at wavelengths A <J 100 /im, and fit by theoretical models 
of graphite and silicate dust in the form of small spheres 
with a distribution of radii (Draine & Lee 1984; Draine 
& Shapiro 1984). Because the opacity of this 'standard' 
dust falls off rather quickly with increasing wavelength for 
A ;> 0.1 /im, attenuated light is significantly reddened. 

The weakness of this method when applied to a new 
situation is the necessary assumption that the same 
extinction-reddening relation holds, even though the dust 
may be of a different species; standard techniques would 



not correct for dust which causes extinction without sig- 
nificant reddening. For the same reason, an effectively uni- 
form intergalactic distribution of non-reddening dust could 
remain undetected by reddening studies such as those of 
Wright k, Malkan (1988) and Cheng, Gaskell & Koratkar 
(1991). 

2. A SPECIFIC MODEL WITH CARBON NEEDLES 

To make this idea concrete, let us begin with the the- 
ory of dust formation. Very small dust grains are believed 
to form in a vapor-solid transition, via a process of crys- 
tal formation; these grains may then coagulate into larger 
ones. In small grain formation, nuclei form first, and then 
they grow as surface nucleation occurs on their faces. Sur- 
face nucleation creates steps which grow along the face, 
adding a new layer and increasing the crystal size. But 
environments of low supersaturation (as may commonly 
occur astrophysically) strongly inhibit surface nucleation. 
In such cases, grains may still grow by the mechanism of 
a 'screw dislocation' (see e.g. Frank 1949, Sears 1955), 
or by growth of a rolled-up platelet (Bacon 1960), form- 
ing one-dimensional 'needles' like those commonly found 
in laboratory experiments (e.g. Nabarro & Jackson 1958.) 

Moreover, needles can grow rapidly where 'spherical' 
dust cannot; thus in some situations needles can out- 
compete spherical dust for available metals. This rea- 
soning led Donn & Sears (1963) to predict that interstel- 
lar dust could be dominated by needle- type dust. These 
predictions were partially borne out by the discovery of 
estatite needles in captured interstellar dust (Bradley, 
Brownlee & Veblen 1983.) The needles were not the 
dominant component, but their discovery does demon- 
strate that vapor-solid transitions occur astrophysically, 
and that astrophysical needles can form by the same mech- 
anism as in laboratories (they contained screw disloca- 
tions.) Conducting needles are physically interesting be- 
cause they act as antennas, absorbing electromagnetic ra- 
diation more effectively than standard dust. Several au- 
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thors have proposed models in which such grains thermal- 
ize the microwave background (see, e.g., Wickramasinghe 
et al. 1975; Wright 1982; Hawkins & Wright 1988). 

I have calculated the extinction cross section for nee- 
dles at A < 10 pm using the 'discrete dipole approxima- 
tion' (see, e.g., Draine 1988) as implemented in the pub- 
licly available DDSCAT package and using the accompa- 
nying graphite dielectric constants e± and eii (see Laor & 
Draine 1993). Following Wickramasinghe & Wallis (1996), 
I assume that the graphite c-axis is perpendicular to the 
needle so that e± applies for an electric field parallel to 
the needle axis (see also Bacon 1960). Figure H shows 
curves for various needle diameters d, 0.02 < d[pm\ < 0.2, 
averaged over incident radiation directions, and averaged 
over an aspect ratio distribution n(L/d) ex (L/g?) -1 with 
4 < L/d < 32. This mass-equipartition distribution was 
chosen to represent a shattering spectrum from longer nee- 
dles; laboratory needles grow up to L/d ~ 1000. The 
maximal L/d is somewhat arbitrary but largely irrelevant 
since the short-wavelength behavior depends only weakly 
on L/d > 8. The results roughly agree with the Mie cal- 
culations of Wickramasinghe and Wallis (1996) which use 
somewhat different optical data. 

Major uncertainties in the opacity model include the 
uncertainties in optical data (see Draine & Lee 1984 for 
discussion), an unknown impurity content in the needles, 
and the unknown needle diameter; 1 the model given is in- 
tended to be suggestive rather than complete. The key 
point is that needles gcncrically have an opacity which is 
higher (kv — 10 5 cm 2 g _1 ) and less wavelength-dependent 
than that of standard dust. 

Several works (Chiao & Wickramasinghe 1972, Ferrara 
et al. 1990 and 1991, Barsella et al. 1989) have studied 
dust ejection from galaxies, all concluding that most spi- 
ral galaxies could eject (spherical) graphite dust. These 
theoretical studies are supported by observations of dust 
well above the gas scale height in galaxies (Ferrara et al. 
1991) and of vertical dust lanes and fingers protruding 
from many spiral galaxies (Sofue, Wakamatsu & Malin 
1994). The high opacity of needles extends over a large 
wavelength range, hence they are strongly affected by ra- 
diation pressure and are even more likely be ejected than 
spherical grains. In a magnetized region, charged grains 
spiral about magnetic field lines. Magnetized gas may es- 
cape from the galaxy via the Parker (1970) instability, or 
grains may escape by diffusing across the lines during the 
fraction of the time that they are uncharged; see, e.g., 
Barsella et al. (1989). Once free of the galaxy 2 , needles 
would rapidly accelerate and could reach distances of 1 
Mpc or more. 

Following Hoyle & Wickramasinghe (1988), we estimate 
the time required for needle ejection and dispersement as 
follows. A grain with length L, cross section d 2 , specific 
gravity p m and opacity k in a anisotropic radiation field 



will attain a terminal velocity v given by equating the 
radiative acceleration kF/c to the deceleration due to vis- 
cous drag of a « (4v 2 /ird)(p ga s / Pm) + 0(u 2 /v 2 ). Here, 
F is the net radiative flux, and u — (3kT gas /mH) and 
pgas are the gas thermal speed and density. Values ap- 
plicable for needles in our Galaxy are k ~ 10 5 cm 2 g _1 , 
p m ~ 2 g_cm- 3 , p g ~ 10~ 24 g cm" 3 , T gas ~ 100K, and 
F/c ~ 10 13 ergs cm 3 . These give a terminal velocity of 
ti~4x 10 5 (c?/0.1/zm) 1 / 2 cm s" 1 and a timescale to escape 
a 100 pc gas layer of - 2.5 x 10 7 (d/0.1^m)- 1 / 2 yr. Out- 
side the gas layer, the needle is subject only to radiation 
pressure. For a rough estimate we assume that the con- 
stant acceleration kF/c ~ 10~ 8 cm s~ 2 acts for the time 
required for the needle to travel a distance equal to the 
galactic diameter. This takes ~ (2Rc/ 'kF) 1 / 2 ~ 8 x 10 7 yr 
for a galaxy of size R ~ 3 x 10 22 cm and leaves the needle 
with velocity v ~ 2.5 x 10 7 cm s — 1 . Such a velocity will 
carry the needle 1 Mpc (twice the mean galaxy separation 
at z = 1) in ~ 4 Gyr. For comparison, the time between 
z = 3 (when dust might be forming) and z = 0.5 (when 
the supernovae are observed) is 5.5 Gyr for f2 = 1 and 
7.3 Gyr for = 0.2. These estimates suggest that radia- 
tion pressure should be able to distribute the dust fairly 
uniformly 3 before z ~ 0.5. 

Dust is known to exist in large quantities (masses ~ 
0.1% of the total galaxy mass are often inferred) in bright, 
high-redshift galaxies (see, e.g., Hughes 1996). These 
galaxies would preferentially eject dust with higher opacity 
at long wavelengths (e.g. needles, or fractal/fluffy grains); 
such grains tend to have a shallower falloff in opacity with 
wavelength, hence redden less than the observed Galac- 
tic dust. This selection affect and the estimation of dust 
escape timescales suggest that if substantial intergalactic 
dust exists, it should be effectively uniform, and redden 
less than standard dust. 

We can compute the optical depth to a given redshift 
due to uniform dust of constant comoving density using 
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Figure [l] shows the integrated optical depth to various red- 
shifts for needles with d = 0.1 pm, for = 0.2, h = 0.65 
and flneedie = 10~ 5 . Using this information, we can calcu- 
late the dust mass necessary to account for the observa- 
tions if ft a = 0. The difference between an ft — 0.2, 17a 
().() model and a model with il = 0.24, ft A = 0.76 (the 
favored fit of R98) is about 0.2 magnitudes at z = 0.7. 
In the d = 0.1 pm needle model this requires Qneedle = 
1.6 x 10~ 5 . Matching an f2 = 1, J7a = universe re- 
quires about 0.5 magnitudes of extinction at z = 0.7 and 

^needle = 4.5 X 10~ 5 . 

A reddening correction based on standard dust proper- 
ties, like that used in R98, would not eliminate this effect. 
R98 effectively estimates extinction using rest-wavelength 

lr The needle diameter is particularly important, but it is difficult justify any a priori estimate of its value. For the sake of the argument at 
hand I take d = 0.1. Note that a distribution of diameters would likely simply be dominated by the low-diameter cutoff. 

2 A grain leaving the disk would be sputtered by gas in the hot galactic halo, but the resulting mass loss is < 20% for the 0.01 /im silicate 
spheres and even less for faster moving or larger grains (Ferrara et. al. 1991), so the effect on the needles would be very small. 

3 Using the model of Hatano et al. (1997), R98 argues that dust confined to galaxies would cause too large a dispersion in supernova fluxes. 
For the dust to create less dispersion than observed, they must merely be 'uniform' enough that a typical line-of-sight (of length ~ cH^ ) 

passes through many clumps of size ~ D and separation ~ A, i.e. that 1/-^/ cHq 1 D 2 / 'A 3 < 1 (the observed dispersion in R98 is 0.21 mag, 
similar to the necessary extinction). This is easily satisfied for needles traveling > 50 kpc from their host galaxies. 
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(after K-correction) B — V color and the Galactic redden- 
ing law. For standard dust this would be reasonable even 
for a cosmological dust distribution, since the reddening 
would still occur the across the redshift-corrected B and 
V frames. But Figure [j] shows that this does not hold 
for needles: the d = 0.1 /im needle distribution only gives 
(B — V) = 0.06Ay up to z — 0.7. The supernova group 
method would if —correct the B and V magnitudes, then 
convert this (rest frame) B — V into an extinction based on 
the Galactic (B — V) — 0.32Ay. It would therefore not be 
surprising for the systematic extinction to go undetected. 

Studies of redshift-dependent reddening (e.g. Wright 
1981, Wright & Malkan 1987, Cheng et al. 1991) in far- 
UV (rest frame) quasar spectra put limits on a uniform 
dust component, but these are most sensitive to high red- 
shifts, at which the needles would not yet have formed and 
uniformly dispersed. In addition, it is clear from Figure 
[l] that for thick needles the flatness of the opacity curve 
would lead to a very small shift in the quasar spectral 
index up to z = 1. 

Another available constraint, the metallicity of Ly-a 
clouds, is probably not relevant; because the dust forma- 
tion and ejection (due to radiation pressure) from galaxies 
is independent of the enrichment mechanism of the clouds 
(presumably population III enrichment or gas 'blowout' 
from galaxies), there is no clear connection between the 
mass of metal gas in the clouds and the mass of needle 
dust in the IGM 4 . 

To estimate the fraction of carbon locked in the nee- 
dle dust, we would like to know Viz at 0.5 <J z <; 3. 
The current value of £lz should be bounded above by the 
metal fraction of large clusters, which are the best avail- 
able approximation to a closed system that is a fair sam- 
ple of the universe. Clusters tend to have ~ 1/2 solar 
metallicity (e.g. Mushotsky et al. 1996), and ~ 10% of 
their mass in gas (e.g. Bludman 1997 for a summary), 
giving Viz <y 10~ 3 . This compares reasonably well with 
an upper bound on universal star density estimated from 
limits on extragalactic starlight (from Peebles 1993) of 
£1* < 0.04: if we extrapolate the Galactic metallicity of 
~ Zq, we find flz ~ ZqVL* <J 4x 10~ 4 . Assuming a current 
Viz ~ 4x 10~ 4 and that metals are created constantly (con- 
servative, given the higher star formation rate at high-z) in 
time from z = 6 we find (for both Vt = 0.2 and Q = l) that 
Vl z {z = 3) - 4 x 10" 5 and fl z {z = 0.5) - 2 x 10" 4 , which 
agrees with recent estimates by Renzini (1998). Even 
such crude approximations are very vulnerable, but sug- 
gest that the needed amount of needle mass is reasonable. 

The needle model is falsifiable in several ways. First, the 
needle opacity spectrum is not perfectly flat, especially for 
small d. Observations over a long wavelength span might 
reveal a redshift-dependent systematic change in certain 
colors. 

Next, the needles take some minimum time to form, then 
more time to achieve a uniform cosmic distribution. Thus 
at high enough redshift the dispersion in supernova bright- 
nesses discussed in R98 appears. Moreover, at z = 1.5 the 
difference between the O = 0.2, fi^ = model with dust 



and the = 0.24, Q\ — 0.76 model sans dust is ~ 0.2 
mag, which should eventually be observable. 

I shall not attempt to address the question of galaxy 
counts here. As commented in R98, grey dust would exac- 
erbate the 'problem' of unexpectedly high galaxies counts 
at high-z, but the magnitude of such an effect would de- 
pend upon the dust density field's redshift evolution, and 
a full discussion of the galaxy count data as a constraint 
on the model (requiring also an understanding of galaxy 
evolution) is beyond the scope of this letter. 

Galactic observations probably cannot disprove the 
model, since needles with the properties most different 
than those of Galactic dust would be ejected with high 
efficiency. Moreover, dust with needle-like characteristics 
may have been detected by COBE (Wright et. al. 1991; 
Reach et al. 1995; Dwek et al. 1997) as a minor 'very cold' 
component of Galactic dust. Such a component is best ex- 
plained by dust with a hitherto-unknown IR emission fea- 
ture, or by fluffy/fractal or needle dust (Wright 1993) and 
could represent a residual needle component with about 
0.02-0.4% of the standard dust mass. 5 

On the other hand, the dust cannot escape from clus- 
ters, which have much higher mass/light ratios, so nee- 
dles formed after the formation of a cluster should re- 
main trapped within. Studies of background quasar counts 
(Bogart & Wagner 1973; Boyle et al. 1988; Romani 
& Maoz 1992), cooling flows (Hu 1992), and IR emis- 
sion (Stickel et al. 1998) of rich clusters indicate extinc- 
tions Ay ~ 0.2 — 0.4 mag and standard dust masses of 
M dust _ 1Q io Denoting by Z cU M cU and M^ s the 

mean cluster metallicity, total mass and gas mass, we can 
estimate the fraction of metals in dust Xcl to be 



M, 



Xcl 



M dust 



cl 



1 cl 



10 x 10~ 5 /0-01 = 0.01, 



using M c i ~ 10 15 Mq. Comparing this to the Xgai < 1 
typical of our Galaxy would indicate dust destruction 
efficiency of <; 99% in clusters. An earlier calcula- 
tion gave Vlneedie/Qz <; 0.1 for the intergalactic needles. 
Assuming the calculated dust destruction, this predicts 

M needle l M dust ^ Q L The needles are about five times a£j 

opaque in optical as standard dust, so this gives an opti- 
cal opacity ratio of ~ 0.5. If these estimates are accurate, 
comparison of nearby cluster supernovae to nearby non- 
cluster supernovae at fixed distance should reveal a mean 
systematic difference of Ay J> 0.03 — 0.06 in fluxes after 
correction for reddening. 6 The Mt. Stromlo Abell cluster 
supernova search (Reiss et al. 1998), currently underway, 
should make such an analysis possible. Note that uncer- 
tainties in the needle opacity relative to standard dust will 
not affect the cluster prediction which (modulo the quan- 
titative uncertainties) should hold unless clusters destroy 
needles more efficiently than standard dust. 

3. CONCLUSIONS 



I have argued that the reduction of supernova fluxes at 
high redshift could be caused by a uniform distribution of 
intergalactic dust. 

4 Of course, if the needles were also assumed to form in the population III objects their density should then relate to the Ly-a metallicity. 

5 The needles absorb about 5x as effectively (per unit mass) in the optical where most Galactic radiation resides, and the 'very cold' 
component emits between .1% and 2% of the total FIR dust luminosity (Reach et al. 1995). 

6 Similar arguments might apply to elliptical galaxies, from which dust ejection is less efficient than from spirals. 
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Both theoretical arguments and observational evidence 
strongly suggest that some dust should be ejected from 
galaxies. Dust with high opacity (especially at long wave- 
lengths where most of the luminosity of high-redshift star- 
burst galaxies resides) would be preferentially ejected. But 
this is exactly the sort of dust which would both redden 
less than standard dust, and require less dust mass to pro- 
duce the observed effect. This letter develops a specific 
model of intergalactic dust composed of carbon needles- 
a theoretically reasonable, even expected, form of carbon 
dust-with conservative properties. The supernova data 
can be explained by a quantity of carbon needles which is 
plausible in light of rough estimates of the universal metal 
abundance. 

Because the dust distribution is effectively uniform, it 
does not induce a dispersion in the supernova magnitudes, 
and because it absorbs more efficiently than standard dust, 
it does not require an unreasonable mass. Finally, because 
the dust is created and ejected by high-z galaxies, it does 
not overly obscure very high redshift galaxies or quasars. 
Thus the key arguments given in R98 against 'grey' dust 
do not apply. The dust of the proposed model should, 
however, provide independent signatures of its existence; 
one is a systematic difference in fluxes between cluster and 
non-cluster supernovae which may be detectable in ongo- 
ing surveys. Finally, the needle model is only one specific 
model for intergalactic dust. Other possible 'dust' types 
are fractal dust (e.g. Wright 1987), platelets (e.g. Donn & 
Sears 1963; Bradley et al. 1983), hollow spheres (Layzer 
& Hively 1973), or hydrogen snowflakes. 

The explanation of reduced supernova flux at high red- 



shift described in this letter depends upon the plausi- 
ble but still speculative assumption that the intergalactic 
dust distribution has significant mass, and is dominated 
by grains with properties exemplified by those of carbon 
needles. The probability that this is the case should be 
weighed against the severity of the demand that the ex- 
planation favored by Riess et al. and Perlmutter et al. 
places on a solution of the vacuum energy (or cosmologi- 
cal constant) problem: the expected value of the vacuum 
energy density at the end of the GUT era must be reduced 
by some as yet unknown process, not to zero, but to a value 
exactly one hundred orders of magnitude smaller. 

This papers has benefited significantly from the com- 
mentary of an anonymous referee. I thank David Layzer, 
George Field, Alyssa Goodman, Bob Kirshner, Ned 
Wright and Saurabh Jha for useful discussions. 

4. NOTE ADDED IN PROOF 

Two recent papers bear upon this Letter. Kochanek 
et al. (1998, |astro-ph/981111l| ) find that dust in early- 
type high z galaxies reddens more than 'standard' dust. 
Needles may escape ellipticals, but the lensing technique 
applied to clusters would be an excellent test of the needle 
model. 

Perlmutter et al. (1998, Ap. J., accepted) find no sta- 
tistically significant different in mean reddening between 
low-z and high-z samples; however, the B — V ~ 0.01 mag 
that my fiducial model predicts fits easily within their l-er 
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Fig. 1. — Top: Carbon needle opacity curves for aspect ratio distribution n oc (L/d)~ 1 and diameters d = 0.02,0.04,0.1,0.2 (im (dotted, 
dashed, solid and triplc-dot-dashcd, respectively). Also included are central U, B, V and I frame wavelengths and a 'Drainc & Lee' dust curve 
(long-dashed) calculated using the Laor & Draine (1993) method. Bottom: Integrated extinction to various redshifts for d = 0.1 fim needles 
for H = 0.2 and Cl neeMe = lO" 5 . 



